Relation between quantum dwell times and flux-flux correlations 
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We examine the connection between the dwell time of a quantum particle in a region of space 
and flux-flux correlations at the boundaries, ft is shown that the first and second moments of a 
flux-flux correlation function which generalizes a previous proposal by Pollak and Miller [E. Pollak 
and W. H. Miller, Phys. Rev. Lett. 53, 115 (1984)], agree with the corresponding moments of the 
dwell-time distribution, whereas the third and higher moments do not. We also discuss operational 
approaches and approximations to measure the flux-flux correlation function and thus the second 
moment of the dwell time, which is shown to be characteristically quantum, and larger than the 
corresponding classical moment even for freely moving particles. 

PACS numbers: 03.65.Xp, 03.65.Ta, 06.30.Ft 



I. INTRODUCTION 

Time observables, i.e., times measured as random vari- 
ables after a microscopic system is prepared, e.g. accord- 
ing to a given wavefunction, are very common in the labo- 
ratories. Examples arc the lifetimes of excited species, or 
the arrival times of ions, atoms or molecules at a scintilla- 
tion detector, a microchannel plate or a laser-illuminated 
region. Nevertheless, incorporating time observables into 
quantum mechanics is a problematic task and it has 
even been discouraged by many physicists influenced by 
Pauli's theorem and the extended notion that "time is 
only a parameter" within the quantum realm. Time ob- 
servables are characterized experimentally, for a given 
preparation of the initial state, by a distribution or by 
its statistical moments, and useful information may be 
extracted from them. The spectacular progress in quan- 
tum state manipulation with laser and magnetic cooling 
techniques emphasizes the need to treat atomic motion 
quantally rather than classically 0, E|, and therefore a 
quantum approach to time quantities. Much work has 
been carried out in the last two decades trying to formu- 
late time observables theoretically and also to connect the 
abstract proposals with actual or idealized experiments. 
Central to these investigations have been the tunneling 
time, the arrival time, and the dwell time 

The dwell time of a particle in a region of space and its 
close relative, the delay time Q, are in particular rather 
fundamental quantities that characterize the duration of 
collision processes, the lifetime of unstable systems, the 
response to perturbations, or the properties of chaotic 
scattering. In addition, the importance of dwell and de- 
lay times is underlined by their relation to the density of 
states, and to the virial expansion in statistical mechan- 
ics @. 

The theory of the quantum dwell time is quite peculiar 
and subtle in several respects To begin with, unlike 
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other time quantities, there has been a broad consensus 
on its operator representation d, Q , 



Tr 



(1) 



where Xfi(t) is the (Heisenberg) projector onto a region 
of space, which we shall limit here to one dimension for 
simplicity, 1Z = {x : x\ < x < X2} , 



Xn{t) = e lRt ' h dx\x)(x\e- im /\ 



(2) 



We shall also assume that the Hamiltonian holds a purely 
continuous spectrum with degenerate (delta-normalized) 
scattering eigenfunctions \4>±k) corresponding to incident 
plane waves | ± k), with energy E = k 2 h 2 /(2m). The 
concept of an average dwell time 7~d(A;) for a finite space 
region in the stationary regime is due to Biittiker 

1 



spac 

m 



r D (k) 



dx\Mx)\, 



(3) 



where j(k) is the incoming flux associated with \4>k)- 

The operator Tjj is positive definite and essentially self- 
adjoint. Moreover, being a "time duration" rather than 
a time instant, To commutes with the Hamiltonian with- 
out conflict with Pauli's theorem, and therefore it can be 
diagonalized in the eigenspace of H. This simplifies the 
derivation of the corresponding quantum dwell-time dis- 
tribution, which, for a state = \ip(t = 0)) is formally 
given by 



n(r)=(^(T D -r)|^). 



(4) 



Following the same manipulation done for the S operator 
in one-dimensional scattering theory, it is convenient to 
define an on-the-energy-shell 2x2 dwell-time matrix T, 
by factoring out an energy delta, 



h\T D \<t>' k )=5(E k -E k ,)^T kk , 



where 



= 6(k - k') and 



Tfefc' = {4>k\Xn\4>k') 



hm 



Ei. — Ei, 



(5) 



(6) 
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In particular, Tu(k) = Tkk- Diagonalization of T for free 
motion leads generically to an interesting quantum pe- 
culiarity (26[: the existence of two different dwell-time 
eigenvalues t±(k) — mL[l ± sm(kL)/kL]/Hk on the en- 
ergy shell, L = x 2 — xi, at variance with the classical 
case. As a consequence, a broader spread (variance) in 
the quantum distribution than the one expected for a 
corresponding ensemble of classical particles will occur. 

In spite of the nice properties of Tq, a direct and suf- 
ficiently non-invasive measurement of the dwell time of a 
quantum particle in a region of space, so that the statis- 
tical moments are produced by averaging over individual 
dwell-time values, is yet to be discovered. If the particle 
is detected (and thus localized) at the entrance of the 
region of interest, its wavefunction is severely modified 
("collapsed"), so that the times elapsed until a further 
detection when it leaves the region do not reproduce the 
ideal dwell-time operator distribution, and depend on the 
details of the localization method. Proposals for oper- 
ational, i.e. measurement-based approaches to traversal 
times based on model detectors studying the effect of lo- 
calization have been discussed by Palao et al. [ll[ and 
by Ruschhaupt jT2j. All operational approaches to the 
quantum dwell time known so far provide in fact, and 
only indirectly, just the average, by deducing it from its 
theoretical relation with some other observable with mea- 
surable average. It is obtained for example by a "Larmor 
clock" , using a weak homogeneous magnetic field in the 
region 1Z and the amount of spin rotations of an incident 
spin-i particle [13, [H, G3 ■ An optical analogue is pro- 
vided by the "Rabi clock" [l5[. It can also be obtained 
by measuring the total absorption if a weak complex ab- 
sorbing potential acts in the region [r| [13, EH- This 
setup could be implemented with cold atoms and lasers 
as described in 

The expression for the average of the dwell time for 
time-dependent scattering processes in terms of the prob- 
ability densities corresponds (in sharp contrast to its sec- 
ond moment, as shown below) to the classical expression 
for an ensemble of classical particles and it reads (20I . |2~0 ] 

/OO [>X2 
dt dx\^j(x,t)\ 2 (7) 
-OO J X\ 

= / dk\(k\^ n )\ 2 r D (k), (8) 
Jo 

where ip(x,t) = J" °° dk (<fik\4>} exp(— ihk 2 t/2m)(j)k{x) is 
the time-dependent wave packet and we assume, here 
and in the rest of the paper, incident wavepackets with 
positive momentum components. To write Eq. ([5]) 
use has been made of the standard scattering relation 
(<f> k \il>) = (k\ip m ), where (x\k) = (2tt)~ 1 ^ exp(ikx), and 
ip m is the freely-moving asymptotic incoming state of if). 
Integrals of the form {7} had been used to define time 
delays by comparing the free motion to that with a scat- 
tering center and taking the limit of infinite volume [22T | . 

For a sample of further theoretical studies on the quan- 
tum dwell time see @, H S3, SI SI, H SI, SI, HE HI 



|32j. A recurrent topic has been its role and decomposi- 
tion in tunneling collisions. Instead, we shall focus here 
on a different, so far overlooked, but rather fundamental 
aspect, namely, the measurability and physical implica- 
tions of its second moment. 

We shall generalize the approach by Pollak and Miller 
[HI], who showed that the average stationary dwell time 
agrees with the first moment of a microcanonical flux-flux 
correlation function (ffcf). We demonstrate that this re- 
lation holds also for the second moment, and we extend 
their analysis to the time-dependent (wavepacket) case. 
The relation fails for third and higher moments and thus 
the ffcf contains only part of the information of the dwell- 
time distribution, although it is certainly the most rele- 
vant. We shall also discuss a possible scheme to measure 
ffef's, thus paving the road towards experimental access 
to quantum features of the dwell time distribution. 



II. STATIONARY FLUX-FLUX CORRELATION 
FUNCTION 

A connection between the average stationary dwell 
time and the first moment of a ffcf has been shown by 
Pollak and Miller [33[. They define a quantum micro- 
canonical ffcf CpM(T,k) = Tr{ Re Cpm (t, k)} by means 
of the operator 

C pm (t, k) = 2nh[J(x2,T)J(x 1 ,0) + J{ Xi ,t)J(x 2 ,0) 
- Jixu^JixuQ) ~ J(x 2 ,T)J(x 2 ,0)]5(E k ~H), (9) 

where J(x,t) = e lHt/h ^[p5(x - x) + S(x - x)p\e- lHt / h is 
the quantum mechanical flux operator in the Heisenberg 
picture, and p and x are the momentum and position 
operators. 

The motivation for this definition stems from classical 
mechanics and can be understood intuitively: Eq. ^ 
counts flux correlations of particles entering 1Z through 
x\ (x 2 ) and leaving it through x 2 (xi) a time r later. 
Moreover, particles may be reflected and may leave the 
region 1Z through its entrance point. This is described by 
the last two terms, where the minus sign compensates for 
the change of sign of a back-moving flux. Note that these 
negative terms lead to a self-correlation contribution that 
diverges for r — > 0. 

We review in the following the derivation of the average 
correlation time and show afterwards that a similar rela- 
tion for the second moment holds. As in the rest of the 
paper, we shall only consider positive incident momenta, 
so that we shall actually deal with Cp~ M , substituting 
S(E k - H) by 5+{E k - H) := 8{E k - H)A + , where A+ is 
the projector onto the subspace of eigenstates of H with 
positive momentum incidence. First of all, we note that, 
by means of the continuity equation, 

4~J(x,t) = ^-J(x,t), (10) 



3 



where p(x, t) = e lHt / h 5(x — x)e lHt / h is the (Heisenberg) 
density operator, Cp M (r, k) can be written as 



was not observed in Ref. [32 
way as above, 



Proceeding in a similar 



C+ M (T,k)=-2nh 



dr 



Xn{r) 



dt 



Xn(t) 



(11) 



By a partial integration and using the Heisenberg equa- 
tion of motion the first moment of the Pollak-Miller cor- 
relation function is given by 



Tr 



dr tC£ m (t, k) 



TV { 2 nil [ drxn(r)^[xn(0), H]S+(E k - H) 



(12) 



Boundary terms of the form \\u\ T —, 00 T 1 X'Ti{ T )il = 1,2, 
are omitted here and in the following. This omission 
can be justified by recalling that physical states must be 
time-dependent and square-integrable so that the contri- 
bution of these terms should vanish when an integration 
over stationary wavefunctions is performed to account 
for the wavepacket dynamics. (In the next section we 
shall discuss explicitly a time dependent version of the 
correlation function.) For potential scattering the prob- 
ability density decays generically as r~ 3 , which assures a 
finite dwell time average, but for free motion it decays as 
t -1 [H]], making td infinite, unless the momentum wave 
function vanishes at k — sufficiently fast as k tends to 
zero 



2G 



Writing the commutator explicitly and using the cyclic 
property of the trace gives 



S+(E k -H). Tr^J diVC+,(T,fc) 

= Tr UttH ( dTXn(r)xn(0)S+(E k -H)\. (15) 



With Tr{- • • } = J dk'{<f>k> | • • • \4>k') and taking "the real 
part of" Cp M , i.e. [Cp~ M + (Cp M Y}/2, we obtain 



Tr Re / drr 2 C+ M (r,fc) = 





Re 


[C 






'xi 


(T 2 


)kk- 



dx \4>k{x)[ 



47r 2 m 2 



dec cj)* k {x)<j)- k (x) 



(16) 



Interestingly, one finds exactly the second moment of T. 
This shows that the relation between dwell times and 
flux-flux correlation functions goes beyond average values 
and that Cp M (r,k) includes the quantum features of a 
dwell time: note that the first summand in Eq. (|16p 
is nothing but (T kk ) 2 , whereas the second summand is 
positive, which allows for a non-zero on-the-energy shell 
dwell-time variance. We insist that the stationary state 
considered has positive momentum, (f>j~(x), k > 0, but 
this second term implies the degenerate partner <p_ k (x) 
as well and is generically non-zero. 

The question arises if these connections hold for the 
other moments of Cp~ M (r, k) . The answer is no, as we will 
show in the next section with a more general approach. 



Tr / drrC+ M (r,k) 



Tr { 2nh dr xr(t) ) xn(0)S + (E k - H) 



(13) 



and integration over r yields the final result of Pollak and 
Miller, 



Tr|^ drrC+ M (r,k) 
27rhTr{xn{0)S + {E k -H)} 



= T 



kk- 



(14) 



Expressing the trace in the basis \(f> k ) gives back the sta- 
tionary dwell time of Eq. i.e. the diagonal element of 
the on-the-energy-shell dwell time operator, T kk . 

Note that in our derivation it becomes obvious that 
the average correlation time is real, in spite of the fact 
that Cp M (r, k) is not selfadjoint. (The discussion of the 
imaginary part of the ffcf in 33] is based on a modified 
version of Eq. ([9]).) 

Next, we will show that the second moment of the 
Pollak-Miller ffcf equals the second moment of T. This 



III. TIME-DEPENDENT FLUX-FLUX 
CORRELATION FUNCTION 

In the following we present a time-dependent version 
of the above flux-flux correlation function and show its 
relation to dwell times. So far, ffcf have been mostly 
considered in chemical physics to define reaction rates 
for microcanonical or canonical ensembles [351 ] - However, 
a physically intuitive time-dependent version can be de- 
fined in terms of the operator 



C(t)= / dt[j(x 2 ,t+T)J(x 1 ,t) + J(x 1 ,t+T)J(x 2 ,t) 
J —oo 

- J( Xl ,t + T)J(xi,t) - J(x 2 ,t + r)J(x 2 ,t)], (17) 



which leads to the flux-flux correlation function 



C(t) = (ReC(r))^, 



(18) 



where the real part is taken to symmetrize the non- 
selfadjoint operator C(t) as before. 

As in the stationary case, Eq. (fig)) counts flux correla- 
tions of particles entering 1Z through x\ or x 2 at a time 
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t and leaving it either through x\ or X2 a time r later. 
Moreover one has to integrate over the entrance time t. 
It is easy to show that the first moment of the classical 
version of Eq. (|17p. where J is replaced by the classical 
dynamical variable of the flux, gives the average of the 
classical dwell time. 

As in Eq. (|lip . we may rewrite C(r) in the form 



C(r) 



dt-^-Xn{x,t + T)^-Xn(x,t). (19) 
dr at 



First of all we note that the ffcf C(r) is not normalized, in 
fact its negative contributions exactly cancel the positive 
ones, 



drC(r) 



dtxn(t)-^Xn(t) = 0- (20) 



Next we derive the average of the time-dependent cor- 
relation function. With a partial integration one finds 



&ttC{t) 



oo poo 



dt / dTXn{x,t+T)—xn{x,t). 

oo JO flr 



(21) 



A second partial integration with respect to t, replacing 
d/dt by d/dr and integrating over r gives 



drrC(r)=/ dtx\{x,t) = dt x n (x, t) = T D , 

J — oo J —oo 

(22) 

where the projector property of xk has been used. 
Eq. ([22| generalizes the result of Pollak and Miller to 
time-dependent dwell times. 

A similar calculation can be performed for the second 
moment of C(r). After three partial integrations with 
vanishing boundary contributions to get rid of the factor 
t 2 one obtains 



/oo />oo 
dt / dTXn{x,t+r)xn(x,t). 
-oo JO 

(23) 



We evaluate the real part of this expression. Making the 
substitutions t + r — > t and r — > — r in the complex 
conjugated term, we find 



Rc 



Att 2 C(t)=TI. 



(24) 



IV. EXAMPLE: FREE MOTION 

In this section we study the simple case of free motion 
to show the relations between dwell times and ffcf. For 
a stationary flux of particles with energy E^, k > 0, de- 
scribed by (j>k{x) = (x\k) = (27r) -1 / 2 e lfcx , the first three 



moments of the ideal dwell-time distribution on the en- 
ergy shell are given by 



T/cfc — 



mL 



(25) 



(T 2 )fcfc 



(T 3 )/c/c 



rn 



2 L 2 / sin 2 (fcL) \ 

(26) 

m 3 L 3 / o sin 2 (fcL)\ 



As proved above, the first two moments agree with the 
corresponding moments of the Pollak-Miller ffcf, but for 
the third moment we obtain 



Tr 



Re / drr 3 C+ M (T,k) 



3r3 r 



m 6 L 
h 3 k 3 



3[l + cos 2 (fcL)] 3sin(2fcL) 



k 2 L 2 



L 3 k 3 



(28) 



In Fig. [T] the first three moments are compared. The 
agreement between (T 3 )/^ and Eq. ([28")) is very good for 
large values of k, but they clearly differ for small k. 

However, the agreement of the first two moments sug- 
gests a similar behavior of II(r) and C(r). To derive 
n(r), we use the fact that the eigenvalues of T are 
t±{k) — mL[l ± sin(kL)/kL]/hk and the corresponding 
eigenstates are \t±(k)) = (\k)±exp[ik(xi+X2)]\ — fc))/\/2- 
Then, for a wavefunction ip(k) := (k\ijj) with only posi- 
tive momentum components one obtains from Eq. Q 



n D (r) 



I V V — 

9 \pi 



3 7=± 



(29) 



where ^J(t) are the solutions of the equation Fj(k) = 
t-fik) — t = 0. In particular, we use the following wave- 



function for the calculation [361 ] . 

%{k) = N(l - e- Qfe2 ) e-( fe ~ fc °) 2 /[ 4 ( Afc ) 2 ]e- lfe;ro e(/c), (30) 

where N is the normalization constant and O(fc) the step 
function. For the free flux-flux correlation function we 
write 

/•oo poo 

C(r)=Re/ dk dk'ip*(k)ip(k')(k\d(T)\k'), (31) 



JO JO 

and Cfcfc'(r) = {k\C{r)\k') in the free case is given by 

d 2 



Ckk> (t ) 



2nhk 
f(HkT/m + L 



S(k - k')—r 2f(hkT/m) - f(tikT/m - L) 
dr z 



(32) 



where 

f(x) = -2e mx2 l {2hT ^ 



2m 



1/2 



+ iTO erfi N^x 



(33) 
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The result is shown in Fig. [2j The ffcf shows a hump 
around the mean dwell time, as expected (the area under 
this hump is 0.9993), but it strongly oscillates for small 
r and diverges for r — * 0. As discussed above, this is due 
to the self-correlation contribution of wavepackets which 
are at X\ or X2 at the times t and t + t without changing 
the direction of motion in between. A similar feature has 
been observed in a traversal-time distribution derived by 
means of a path integral approach [37l |. 

In contrast, II (r) behaves regularly for t — > 0, but 
shows peaks in the region of the hump. This is because 
the denominator of Eq. (j29|) becomes zero if the slope 
of the eigenvalues i±(fe) is zero, which occurs at every 
crossing point of t + {k) and t_(fc). 

Another dwell time distribution for free motion can 
be defined based on the operator to = n%L/\p\, obtained 
heuristically from quantization of the classical dwell time. 
The eigenf unctions of this operator are momentum eigen- 
functions, |±fc), fc > 0, and the corresponding eigenvalues 
are twofold degenerate and equal to the classical time, 
mL/h\k\. The distribution of dwell times for this oper- 
ator, as always for positive-momentum states, is given 
by 



mL 



(mL\ 



(34) 



The distribution 7t(t) agrees with C(r) in the region 
of the dwell-time peak and tends to zero for r — * 0. 
However, it does not show the resonance peaks of II(r). 
The on-the-energy-shell version of to, t, is also worth 
examining. By factoring out an energy delta function 
as in Eq. (JSJ) we get for a plane wave |fc) the average 
tfefc = mL/(hk), which is equal to T^, but the second 
moment differs, {t 2 ) kk = (t fefc ) 2 = (T fefc ) 2 < (T 2 ) fefe , see 
Fig. 1; in other words, the variance on the energy shell is 
zero since only one eigenvalue is possible for t. Contrast 
this with the extra term in Eq. (|26p which again empha- 
sizes the non-classicality of the dwell-time operator Tjj 
and its quantum fluctuation. 



V. DISCUSSION 

We have demonstrated that the relationship found by 
Pollak and Miller [33| between the first moment of a dis- 
tribution ffcf and the average stationary dwell time is 
also valid for the second moment and for flux-flux cor- 
relations of wavepackets. On the other hand, this rela- 
tionship is not valid for the third moment. While this 
brings dwell-time information closer to experimental re- 
alization, the difficulty is translated onto the measure- 
ment of the ffcf, which is not necessarily an easy task. 
The simplest approximation is to substitute the expecta- 
tion of the product of two flux operators by the product 
of their expectation values (the product of the current 
densities). Using the wave packet of Eq. (|3U|) . we have 
compared the times obtained with the full expression (Tig)) 
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FIG. 1: Comparison of the first three moments: Tkk, [T 2 )kk 
and (T 3 )^^ (dotted-dashed line) with the corresponding mo- 
ments of the flux-flux correlation function, for a free-motion 
stationary state with fixed k. (Tkk) 2 is also shown (dotted 
line), h = m = 1 and L = 3. 
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FIG. 2: Comparison of dwell time distribution II(r) 
(dashed line) and flux-flux correlation function C(t) (solid 
line) for free motion and L — 50. Furthermore, the 
alternative free-motion dwell-time distribution 7t(t) = 
mL\Tp(rnL/ (hr))\ 2 / (hr 2 ) is plotted (circles). The inset shows 
the momentum distribution according to Eq. (|30[1 with fco = 2, 
Afc = 0.4, a = 0.5 and the eigenvalues t±(k). We set 
H = m = 1 and |a;o| large enough to avoid overlap of the 
initial state with the space region 1Z. 



and with this approximation in Fig. [3] The two results 
approach as — > 0. Other factors that improve the 
approximation are the increase of L and/or of fco- 

This approximation can in fact be improved systemat- 
ically, still making use of ordinary current densities, as 
follows: First we decompose J(xi, t+r)J(xj, t) by means 
of the resolution of the identity 



1 

P 



P + Q, 



(35) 
(36) 
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FIG. 3: Comparison of the relative error of (Tb) using the 
approximation Cq(t) instead of C(r) for free motion, a — 0.5, 
H = m = l and L = 100. 

so that 

J(xi,t + T)J(xj,t) = J{x t ,t + T){\i))(^\ + Q)J(x j ,t). 

(37) 

It is useful to decompose Q further in terms of a basis of 
states orthogonal to \rp) and to each other, )}> 

3 

that could be generated systematically e.g., by means of a 
Gram-Schmidt orthogonalization. Now we can split Eq. 
(fT7|) as follows 

C(r)=C (r)+d 1 (T), (39) 

where 

/oo 
dt [J(x a ,* + r)|^)(V'|J(a! 1 ,i) 
-oo 

+ J( Xl ,t + T)\lP){i}\J(x 2 ,t) 

- J{x U t + T)\iP){^\J{ Xl ,t) 

- J(x 2 ,t + r)\iP)(iP\J(x 2 ,t)], (40) 



Similarly, we define C(r) = Cq(t) + Ci(r) by taking the 
real real part of (^|Co(r) + Ci(t)\i/j). Co is the zeroth 
order approximation discussed before and only involves 
ordinary, measurable current densities The non- 

diagonal terms from C\, (tp\J(xi,t)\tp^}('ip^\J(xj ,t + 
r)\ip) can also be related to ordinary fluxes by means 
of the auxiliary states 

IV-i) - + 

m = (42) 

since one easily finds that 

ty\T{x,t)\ii>f) = -(^JM)!^) - ~{ih\J(x,t)\i>2) 

- ^ 3 \J{x,t)\i> 3 ) + % -{^\J{x,t)\^) 

- fa 2 \J(x,t)\i> 2 ). (43) 

To summarize, the present paper provides a route of 
access to the second moment of the quantum dwell time 
through flux-flux correlation functions. This is inter- 
esting because the second moment is characteristically 
quantum and, unlike the first moment, it differs struc- 
turally from (and is larger than) the corresponding clas- 
sical quantity: on the energy shell the dwell time shows a 
quantum fluctuation (non-zero variance) which vanishes 
classically. While the analysis of the quantum dwell time 
has been mostly limited to its average in the existing 
studies, the present results motivate further research on 
the role played by the second moment of the dwell time 
in fields such as lifetime fluctuations, chaotic systems, 
conductivity, or time-frequency metrology. For exam- 
ple, is the second moment affecting the quality factor of 
atomic clocks? The time spent by the atom in a spa- 
tial region determines their stability, which increases in 
principle for slower atoms, but quantum motion effects 
have been shown to become more and more relevant for 
decreasing velocities [3|| |39(. This an other intriguing 
questions on the quantum dwell time are left for sepa- 
rate studies. 



/oo 
dt [j(x 2 ,< + r)|Vf)(Vf \J(xi,t) 
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